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Abstract 

Background: Recent studies have demonstrated that natural transformation and the formation of highly structured 
populations in bacteria are interconnected. In spite of growing evidence about this connection, little is known 
about the dynamics of natural transformation in spatially structured bacterial populations. 

Results: In this work, we model the interdependency between the dynamics of the bacterial gene pool and those 
of environmental DNA in space to dissect the effect of transformation on adaptation. Our model reveals that even 
with only a single locus under consideration, transformation with a free DNA fragment pool results in complex 
adaptation dynamics that do not emerge in previous models focusing only on the gene shuffling effect of 
transformation at multiple loci. We demonstrate how spatial restriction on population growth and DNA diffusion 
in the environment affect the impact of transformation on adaptation. We found that in structured bacterial 
populations intermediate DNA diffusion rates predominantly cause transformation to impede adaptation by 
spreading deleterious alleles in the population. 

Conclusion: Overall, our model highlights distinctive evolutionary consequences of bacterial transformation in 
spatially restricted compared to planktonic bacterial populations. 
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Background 

Bacterial populations often live in highly structured and 
dense populations, known as biofilms, in nature [1]. Bio- 
films are stable assemblages of bacterial cells that stick 
to a surface and are embedded in a self-secreted matrix 
[1-3]. The proximity between cells is advantageous to 
the bacteria as it allows them to respond to environmen- 
tal stresses more effectively, coordinate their behavior 
better and gain access to nutrients and public goods. 
The proximity can, however, be disadvantageous as it 
intensifies the competition for nutrients and other re- 
sources between the bacteria [2,4]. 

In natural habitats, formation of structured bacterial 
communities can occur concurrently with colonization 
of new environments and bacterial adaptation. In the 
course of colonization, bacteria may be under natural 
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selection to adapt to their new environments. One way 
in which the necessary genetic variation for natural se- 
lection to act upon can be generated is natural trans- 
formation, the active uptake of naked DNA fragments 
from the environment and their subsequent integration 
into the genome. Natural transformation has been re- 
ported in many species of bacteria, including some im- 
portant pathogens like Streptococcus pneumoniae and 
Helicobacter pylori [5] . 

There is growing evidence that biofilm formation and 
natural transformation are interconnected processes 
reviewed in [6]. The physical proximity of bacteria in 
space seems to facilitate horizontal gene transfer be- 
tween them, as shown in 5. pneumoniae [7] and H. pylori 
[8] . Further, protein factors that induce natural transform- 
ation were found to also trigger biofilm formation, leading 
to an increased rate of transformation in the biofilm com- 
pared to the planktonic state in naturally transforming 



© 2014 Moradigaravand and Engelstadter; licensee BioMed Central ttd. This is an Open Access article distributed under the 
terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by/4.0), which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons 
Public Domain Dedication waiver (http://creativecommons.0rg/publicdomain/zero/l.O/) applies to the data made available in 
this article, unless othen/vise stated. 



Moradigaravand and Engelstadter BMC Evolutionary Biology 2014, 14:141 
http://www.biomedcentral.com/1471-2148/14/141 



Page 2 of 9 



bacteria [9-11]. Finally, extracellular DNA is known to be 
abundant in biofilms and also to be involved in biofilm 
formation [12,13]. 

Transformation has the potential to facilitate adapta- 
tion, owing to the creation of co-adapted combinations 
of genes [14-18] but also the acquisition of beneficial 
genes. However, some hypotheses for the evolution of 
natural transformation emphasize physiological advan- 
tages to transformation and regard any genetic effect as 
the by-product of these physiological benefits [19,20]. 
Natural transformation can also become costly depend- 
ing on the composition of DNA fragments in the envir- 
onment [15,21]: while acquiring new genes can enhance 
the fitness of the adapting bacteria, transformation with 
defective DNA that originate from bacteria that are 
selected against by natural selection can impede adapta- 
tion. The interplay between costs and benefits of trans- 
formation and bacterial adaptation has been investigated 
in some previous studies [15,21,22]. These studies, how- 
ever, assumed that bacteria live in a planktonic state, and 
therefore, neglected the effect of spatial structure of bac- 
terial populations. Since spatial restriction on population 
growth limits the availability of free DNA to bacteria, we 
expect the impact of transformation on population dy- 
namics to be different in spatially restricted compared to 
the planktonic populations. 

To study the connection between transformation with 
external DNA and structured population growth in bac- 
teria, we developed a mathematical model for adaptation 
of a bacterial population during growth across space. 
The bacterial population was assumed to undergo trans- 
formation with DNA fragments released by dead bac- 
teria. This eco-evolutionary model incorporates several 
ecological and population genetic factors, such as bacter- 
ial migration and DNA diffusion, that govern the rate of 



both population expansion and adaptation in space. We 
demonstrate how the impact of transformation on adap- 
tation varies in different parts of a growing population, 
depending on DNA diffusion, bacterial migration and 
the source of free DNA in the environment. 

Results 

The dynamics of bacterial growth and DNA diffusion 

We first investigated the dynamics of a bacterial popula- 
tion expanding in space in the absence of adaptation. 
Figure 1 shows how a single bacterium gives rise to a 
bacterial colony across space over time. The population 
expands symmetrically because bacteria invade adjacent 
grid cells at the same rate (Figure 1A,B and C). Popula- 
tion growth stops when all grid cells are filled (Figure IC). 
Depending on the migration and growth rates, bacteria 
may migrate to adjacent cells before the carrying capacity 
of a focal grid cell is reached. In this case, a less compact 
population structure is formed over time (not shown). Not 
surprisingly, high levels of bacterial migration lead to 
dynamics resembling those in a well-mbced planlctonic 
population. 

Upon death, bacteria release DNA fragments into the 
environment. In the absence of DNA diffusion, the fi"ee 
DNA remains at the release location and accumulates 
with time. Since the free DNA has more time to amass 
in the initial grid cell, the DNA concentration is higher 
there (Figure 1D,E and F). In our simulations, there is 
always some levels of DNA diffusion and/or DNA deg- 
radation so that free DNA does not indefinitely accumu- 
late in the environment. 

Adaptation with constant selection 

We next considered a growing bacterial population that 
adapts to a new environment under a constant selection 
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Figure 1 Spatial distribution of a growing bacterial population and DNA fragment pool in a 21X21 grid space over time at three time 
points: A,D) 7=200, B,E) 7 = 400 and C,F) 7=500. Simulations start witli a bacterium of genotype 0 located at tlie center of the space. 
Interpolation of order 1 was used to generate continuous plots from discrete data. Parameters take values /• = 0.1, K = 100, m = 10" ^ i/ = 10"^ 
5 = 0, u = 0, /J = 0, g = 0, c = 0. 
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regime where allele 1 has a constant selective advantage 
s over allele 0. The population starts from a bacterium 
of genotype 0, which gives rise to an expanding popula- 
tion. As the population grows, mutation creates bacteria 
of genotype 1, which in the absence of transformation 
go quickly to fixation (Figure 2A). Transformation delays 
the fixation of the beneficial allele in the population and 
impedes adaptation (Figure 2B). This occurs because the 
DNA pool during adaptation is overrepresented by the 
free DNA of type 0 released from the less fit bacteria of 
genotype 0. This so-called 'bad genes effect' of transform- 
ation, also described in prior studies [15,21], is more pro- 
nounced in central layers of the population. This leads to 
the formation of a frequency gradient of bacteria of geno- 
type 1 across different layers of the growing population 
(Figure 2B). In this gradient, the frequency of bacteria of 
genotype 1 is lower in the central layers of the population 
in the short term, because of the higher concentration of 
the free DNA of type 0. 

In order to better understand the impact of migration 
on the evolutionary dynamics of our model, it may be 
helpful to first consider the two extreme ends of the 
spectrum of migration rates. In the complete absence of 
bacterial migration, there is only growth and adaptation 
in the initial grid cell, and this single subpopulation be- 
haves as a planktonic population. Conversely, when mi- 
gration rates are extremely high, the entire population 
becomes well mixed and allele frequencies will be the 
same in all grid cells. Both of these extreme cases can 
thus be understood from previous studies examining the 
impact of natural transformation in planktonic popula- 
tions where the "bad genes" effect was previously de- 
scribed [15,21]. 

To understand how population structure affects the 
"bad genes" effect, we measured the overall frequency of 
bacteria of genotype 1 across the whole population at 
one specific time point, for different intermediate migra- 
tion rates as well as different DNA diffusion rates and 



transformation rates (Figure 3). As anticipated, the bad 
genes effect becomes more pronounced with increasing 
transformation rate. With decreasing bacterial migration 
rate, the "bad genes" effect is mitigated and adaptation 
proceeds faster. This is because the lower the migration 
rate, the longer it takes until new grid cells are colo- 
nized. This delayed colonization of new grid cells means 
that migrants spreading in new grid cells will be better 
adapted to start with and also encounter fewer deleteri- 
ous alleles in the free DNA pool. With low migration 
rates adaptation is therefore impeded by the "bad genes" 
effect in grid cells close to the initial one but proceeds at 
an increasingly faster pace at the frontier of the spatially 
expanding population. Only at extremely low migration 
rates (results not shown in Figure 3) is adaptation again 
impeded because the bacteria fail to colonize any grid 
cell other than the initial one over the time period under 
consideration. 

DNA diffusion potentially has two opposing effects on 
the effect of transformation on adaptation. On the one 
hand, DNA diffusion spreads the free DNA of type 0 in 
the population, and thus, impedes adaptation. On the 
other hand, DNA diffusion depletes free DNA of type 0 
and thus mitigates the bad genes effect. These opposing 
effects can result in a non-monotonic impact of DNA 
diffusion on the transformation impact, particularly with 
low bacterial migration rates (Figure 3B). High levels of 
migration make the non-monotonic effect of DNA diffu- 
sion disappear since in this case, the bacterial population 
behaves like a planktonic population in which the bad 
genes effect is pronounced. At very high DNA diffusion 
rates, the free DNA becomes unavailable to the bacteria, 
and thus the bad genes effect disappears. 

In the Appendix, we present the results of the simula- 
tions in the presence of different levels of DNA degrad- 
ation. Our results demonstrate that DNA degradation 
monotonically reduces the bad genes effect but a rela- 
tively high DNA decay rate (greater than 1/100 of the 
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Figure 2 The number of bacteria of genotype 0 and 1 in different grid cells of the equator section of a 21 x 21 space at one time 
point during adaptation in the absence (A) and presence (B) of transformation. Simulations start with a bacterium of genotype 0 located at 
the center of the space. Red and blue curves correspond to the bacterial numbers of genotype 0 and 1, respectively, and curves corresponding 
to more centrally located grid cells are shown with stronger colors. Parameters take values /• = 0.1, K= 100, m= 10" v= 10""*, s = 0.01, 
(A) (v = Oand (6) (j= 10"^ /j= 10~^ g = 0.1,c = 0. 
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Figure 3 The total frequency of bacteria of genotype 1 in transforming populations growing in a 21 x 21 space at time point 5000 
with different values for bacterial migration and DNA diffusion rates and with two transformation rates: A) u = 10"^ and B) u= 10"^. 

Simulations start with a bacterium of genotype 0 located at the center of the space. Each plot shows the frequencies of allele 1 for 20 parameter 
values of the bacterial migration rate and the DNA diffusion rate. Other parameters take the values r=0.1, K= 100, m = 10"^, 5 = 0.01, ^= 10"^, 
c = 0. 



bacteria growth rate) is required to completely remove the 
effect of transformation (Additional file 1: Figure SI). 
Thus, the dynamics of transformation mentioned above 
are conceivable to be observed in natural circumstances. 

Adaptation with density-dependent selection 

In spatially structured bacterial populations, selection is 
often expected to be heterogeneous across bacteria oc- 
cupying different locations, depending on the density of 
the surrounding subpopulations. For example, restriction 
on the diffusion of an antibiotic makes the interior bac- 
teria less susceptible to the antibiotic, thus decreasing 
the strength of selection for resistance genes. In an ex- 
treme case where no antibiotics reach the interior, there 
might even be selection against resistance genes due to 
costs of resistance, i.e., not only the strength but also the 
direction of selection might depend on the density of 
the surrounding subpopulations. Similar effects may 
arise if diffusion of nutrients is restricted by the bacteria. 



To emulate such effects, we next consider the dynam- 
ics of adapting populations under selection regimes that 
are dependent on population density in the vicinity of a 
focal grid cell. We consider the effect of transformation 
on the distribution of genotype 1 in space in three sce- 
narios. The first scenario is the density-independent re- 
gime where genotype 0 is always disfavored by selection. 
As anticipated in this case, transformation reduces the 
frequency of genotype 1 especially in the center of 
the transforming population during adaptation, due to 
the bad genes effect (Figure 4A). 

The second scenario describes a density dependent re- 
gime where genotype 0 is selected against in an empty 
space but becomes selectively neutral in a full space. In 
this case, genotype 1 becomes more frequent at the less 
dense frontier parts of the population in the absence of 
transformation. DNA diffusion mixes free DNA frag- 
ments of type 0 and type 1 that are over-represented in 
the interior and frontier parts of the population, respec- 
tively. Transformation then reduces the number of 
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Figure 4 The frequency of bacteria of genotype 1 across the equator section of a 21 x 21 grid, for non-transforming (solid curves) and 
transforming (dashed curves) populations growing under different density dependent selection regimes: A) 5^ = 0.01, 5^=0.01, B) Sg = 
0.01, Sf=0 and C) Se = -0.01, Sf=0.01. The curves show the frequencies of allele 1 at time points: A) T= 1000, B) T= 1250 and C) 7=5000. 
(As the speed of adaptation is different in the different selection regimes, different time points have been used to informatively represent the 
dynamics in each regime.) Simulations start from a bacterium of type 0 at the center of the space. Interpolation was used to obtain a continuous 
distribution curve from discrete data points. Other parameters take the values r = 0.1, K= 100, m = 10^^, u= 10~^,/j= 10^^ g= 10~\ c = 0. 
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bacteria of genotype 1 in the frontier of the population 
but increases them in the internal part (Figure 4B). 

The third scenario features a density dependent regime 
where genotype 0 is favored in a free space but selected 
against as population density increases. This results in 
an overrepresentation of bacteria of genotype 1 in the 
interior of the population (Figure 4C). Since genotype 
0 becomes more disfavored as population density in- 
creases, the free DNA of type 0 becomes frequent in the 
environment. Transformation with free DNA of type 0 
lowers the frequency of genotype 1 in the population 
(Figure 4C). Overall, the results of this section indicate 
that in addition to a constant selection regime, trans- 
formation also affects adaptation when the population is 
subject to heterogeneous selective forces. 

Competition between two species 

Thus far, we have assumed that one bacterial species col- 
onizes an empty environment, and therefore, free DNA 
originates from a single colony. However, in natural situ- 
ations bacteria may face already adapted competitor spe- 
cies, so that transforming bacteria are exposed to the 



exogenous DNA fragments originating from the adapted 
species. To examine how transformation affects adapta- 
tion in this case, we study the competition between a 
non-transforming adapted population and an adapting 
population, which potentially undergoes transformation. 
The two populations grow across space under a constant 
selection regime. The non-transforming population car- 
ries allele 1 and is thus assumed to be well adapted to 
the environment; no mutations occur in this population. 
The adapted population therefore supplies the DNA 
fragment pool with only the beneficial allele 1. The 
adapting population consists of two genotypes 0 and 1 
and undergoes mutation and potentially transformation. 
This population has a higher growth rate than the non- 
transforming population. However, bacteria of genotype 
0 in the adapting population are selected against and 
are disfavored in competition with the adapted non- 
transforming bacteria. 

The two populations grow from different parts of the 
grid and compete with each other. Figure 5 reveals that 
in the absence of transformation, the adapted species is 
initially favored and occupies most of the space. This 
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Figure 5 The number of bacteria along the equator layer of a 5 x 100 space in two adapting and adapted bacterial populations at 
three transformation rates for the adapting type: u = 0, u = 10"^ and u= 10"^. The non-transforming population (brown) has already 
adapted to the environment and grows from a bacterium of genotype 1 located at position (3,80). The adapting population consists of two 
genotypes (0: blue, 1 : red) and grows from a bacterium located at position (3,20). The frequencies of different genotypes are taken at time points 
500, 1000, 2000 and 5000. Parameters take value: /■= 0.1 for the adapted population and r= 0.105 for the adapting population, /(= 100, 
m = 10"^, 5 = 0.01, v=10"^(j =10"*^, g= 10^ \c = 0. 
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happens because of selection against genotype 0 in the 
adapting population. As adaptation proceeds, the adapt- 
ing population outcompetes the non-transforming popu- 
lation, owing to its higher growth rate. 

Transformation has two opposing effects on adaptation 
in the adapting population: transformation with the free 
DNA of type 0 released by the adapting species impedes 
adaptation (the bad genes effect), but transformation with 
the free DNA of type 1, which is released by the adapted 
species and diffuses in space, facilitates adaptation. These 
opposing effects result in varying spatiogenetic patterns in 
different parts of the transforming population over time 
(Figure 5). The bacteria at the frontier of the transforming 
population are the first to become exposed to the benefi- 
cial free DNA of type 1, and therefore, undergo adaptation 
faster than the internal bacteria, which are predominantly 
exposed to the conspecific DNA of type 0. Adaptation 
continues in the transforming population and more bac- 
teria of genotype 1 are formed. In the long term, owing to 
their higher growth rate, the transforming population 
gradually out-competes the non-transforming type. 

With a high compared to a low transformation rate, 
although transformation generates more bacteria of 
genotype 1 in the frontline of the population, the bad 
genes effect also becomes more pronounced and this 
outweighs the benefit of transformation. As a result, 
adaptation is considerably decelerated in the transform- 
ing population (Figure 5). Altogether we conclude that 
transformation allows bacteria to take advantage of the 
beneficial free DNA, especially at the frontier of the 
population, and this can outweigh the bad genes effect. 

Discussion 

Lysed bacteria are known to supply the free DNA pool 
in the environment [7,23-25]. This is thought to give rise 
to the "bad genes effect" of natural transformation, ac- 
cording to which there is an overrepresentation of dele- 
terious mutations released by selectively dead cells into 
the DNA pool. Previous models of planktonic bacteria 
considering the dynamics of a free DNA pool have 
shown that this effect results in a fitness cost of natural 
transformation that often overcomes potential benefits 
of transformation [16,21]. Our model shows that spatial 
structure of bacterial populations in combination with 
DNA diffusion strongly affects the bad genes effect. Spe- 
cifically, we found that the bad genes effect is mostly 
mitigated for intermediate bacterial migration and either 
very high or low DNA diffusion rates. 

Our model rests on the assumption that the genetic 
composition of the DNA fragment pool can be different 
from that of the cell population. This is justified because 
free DNA is known to persist in the environment for 
long periods of time and can thus to some extent retain 
a record of genes from previous generations [23]. In a 



spatial context, this also means that DNA diffusion has 
the potential to disseminate genes within the population, 
thereby enabling gene flow without actual bacterial mi- 
gration. However, DNA diffusion may be restricted to 
varying degrees. For instance, extracellular DNA is known 
to serve as an abundant structural polymer in biofilms and 
to be required for the integrity and stability of the extra- 
cellular matrix during biofilm formation [12,13]. Besides, 
bacteria can utilize free DNA as a potential nutrient 
[20,26]. These and other factors can potentially limit the 
diffusion of free DNA in the environment, localize the ef- 
fect of transformation and consequently reduce the extent 
to which transformation increases the level of genotypic 
homogeneity across the population. 

Beyond a transient effect of slowing down adaptation, 
transformation may also lead to a persistent reduction in 
mean fitness of the population when there is spatial het- 
erogeneity, caused by the biofilm structure, in selection 
pressure. One possible scenario that we examined as- 
sumes that selection for an allele becomes weaker or even 
negative at the center of a growing population. A case in 
point is the situation seen in biofilms where bacteria in 
the interior are often immune to the action of antibiotics 
(due to low accessibility of antibiotics [27] or reduced me- 
tabolism caused by nutrients depletion [28,29]), so that re- 
sistance genes are only selected for at the exterior of the 
biofilm. Such heterogeneous selection can entail a migra- 
tion load, i.e. a reduction in the mean fitness of the popu- 
lation due to migration of non-adapted individuals. Our 
results show that transformation can produce an add- 
itional fitness load incurred through incorporating DNA 
that has diffused from parts of the population that are 
under a different selection pressure. 

We also explored a situation where transformation can 
become beneficial during adaptation due to the acquisition 
of free DNA fragments originating from other, already 
adapted populations in the environment. For example, 
consider a situation where a pathogen invades and adapts 
to a complex environment like the nasopharynx where the 
pathogen needs to compete with other bacteria and is ex- 
posed to environmental stresses, such as drug treatment. 
Our results suggest that in a multispecies environment 
transformation becomes advantageous by allowing bac- 
teria to take advantage of the exogenous DNA in the 
shared DNA pool that originated from already adapted 
species. In line with this, prior studies have shown that the 
rate of transformation increases during colonization of the 
nasopharynx by multiple strains of pneumococci [30,31]. 
This effect of free DNA originating from already adapted 
populations resembles the proposed benefit of transform- 
ation in using beneficial genes that the DNA pool retains 
from the past ('the genetic time travel effect') [32]. Simi- 
larly, SzoUosi et al. [33] suggested that horizontal gene 
transfer through bacterial migration may be beneficial by 
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enabling bacteria to acquire genes from other spatially dis- 
tinct populations to readapt faster to changing conditions. 

Our model does not include evolution at multiple loci, 
and therefore, does not involve any potential benefit of 
transformation due to the break up of non-random asso- 
ciations between alleles. Recombination can be beneficial 
during adaptation by bringing together beneficial muta- 
tions that are likely to arise on different backgrounds 
in finite populations (the Fisher-MuUer effect) [34,35]. 
Martens and Hallatschek [36] showed that recombin- 
ation accelerates adaptation in spatially structured popu- 
lations through the Fisher-MuUer effect, but their model 
did not incorporate a free DNA pool. As an additional 
feature of bacterial recombination, our results suggest 
that transformation with the free DNA pool potentially 
affects the fixation of beneficial alleles in structured pop- 
ulations. These findings highlight the importance of con- 
sidering the dynamics of the free DNA pool to fully 
account for the dynamics of recombination in bacteria. 

It should be noted that in our simulations we neglected 
the possibility that genes involved in DNA uptake and 
recombination themselves are subject to homologous 
recombination. In particular, transformation with non- 
functional alleles of such genes released into the environ- 
ment by the non-competent bacteria will produce a self- 
reinforcing drive against transformation so that in the ab- 
sence of any benefits of transformation the competent 
genotype will ultimately become extinct [37] . 

Although our model aims to provide insights into the 
dynamics of spatially structured bacterial populations 
undergoing transformation, our model is not specifically 
tailored to describe bacterial biofilms. Several previous 
models provide a much more reaUstic description of bio- 
film growth. For example, these models explicitly track 
the concentration of different nutrients and growth lim- 
iting/stimulating compounds in different layers of the 
biofilm and account for heterogeneity in position, bio- 
physical and physiological properties of bacteria during 
biofilm formation over time see [38] for a review. In- 
corporating evolutionary dynamics such as those involv- 
ing natural transformation into such individual-based 
biofilm models represents a substantial conceptual and 
computational challenge for future studies. 

Conclusions 

This study was an attempt to elucidate the connection 
between horizontal gene transfer and formation of struc- 
tured population in bacteria. There is growing interest 
among molecular microbiologists to understand this 
connection as researchers have realized that most bac- 
teria live in highly structured populations and horizontal 
gene transfer plays a major role in determining the bio- 
logical repertoire of microorganisms. Experimental tests 
of the above models are needed to gain a better 



understanding of the evolutionary dynamics of bacteria 
during colonization of new environments. 

Methods 

In order to model evolution in a spatially structured popula- 
tion of bacteria, we consider a two-dimensional Lx P grid of 
cells. In each of these grid cells, the bacteria grow at a 
growth rate r that is restricted by the same carrying capacity 
K. The bacteria migrate to adjacent grid cells at a rate v, and 
we consider all bacteria in the 8-cell neighborhood as adja- 
cent. There is no migration beyond the borders of the grid. 

We assume each bacterium has a single locus with 
two alleles (0 and 1) and the bacteria undergo forward 
(1 0) and backward (0 —>■ 1) mutations at the same 
rate ^. The bacteria die at a certain rate that depends on 
their allelic state (see below). Upon death, each bacter- 
ium releases its genomic allele into the medium (initially 
the same subpopulation). In the main text, we assume 
that DNA persists in the environment (no degradation) 
but may diffuse between grid cells at a rate g. In contrast 
to bacterial migration, we assume that the DNA may dif- 
fuse out of the grid of subpopulations. 

The bacteria undergo natural transformation at a rate u 
per bacterium and encounter DNA molecules present in 
the same grid cell. An imported DNA fragment replaces 
the genomic allele through homologous recombination. 
This may lead to the conversion of a genotype to the other 
type (provided the allele that is taken up is different from 
the original one), and it also leads to a reduction of frag- 
ments in the pool of free DNA. Besides DNA diffusion 
and transformation, other environmental factors can limit 
the availability of environmental DNA, such as DNA deg- 
radation and DNA consumption by bacteria as a nutrient 
(see Discussion). To account for this, we include a rate - 
c - at which DNA becomes unavailable to the bacteria. 
The effect of this parameter is studied in the Appendix. 

Taken together, the dynamics of the bacterial popula- 
tion and the DNA pool are described by the following 
system of differential equations: 



x\ = rx\ ( 1- 5^ J -mx\ + V I <-«('./)4 I + 

6 N(,,,) J 

d\ = mx\-ud\(x«+x]\+g[ J2 4?-K -^4 

(1) 
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Table 1 Parameters used in the model 



Parameter Description Unit 

L and P Numbers of grid cells in two dimensional space 

K Carrying capacity of each grid cell - 

r Growth rate Per bacterium per unit time 

^ Mutation rate Per bacterium per unit time 

m Intrinsic death rate Per bacterium per unit time 

u Transformation rate Per bacterium per DNA fragment per unit time 

V Bacterial migration rate Per bacterium per unit time 

g DNA diffusion rate Per DNA fragment per unit time 

s Selection coefficient Per bacterium per unit time 

c DNA decay rate Per DNA fragment per unit time 



In these equations, x and d denote the number of bac- 
teria and DNA fragments, respectively. Subscripts indi- 
cate the grid cell and superscripts indicate the allelic 
state. The term N{i,j) denotes the set of neighboring 
grid cells of the grid cell with the coordinates / and / and 
the term n{i,f) represents the number of these neighbor- 
ing grid cells (i.e., = \N{i,i)\). 

Natural selection operates through differences in death 
rates between the two bacterial genotypes. The total death 
rate is composed of an intrinsic death rate m that is the 
same for all bacteria and a selection coefficient that can 
depend on both genotype and population density. Specif- 
ically, we consider two different selection regimes. In the 
first regime, the allele 0 entails a death rate that is in- 
creased by a constant amount s. Although genotype 0 may 
still grow and establish a population, selection thus re- 
duces the net growth rate of genotype 0 (r-m-s) by an 
amount s compared to the favoured genotype 1. Moreover, 
selection also decreases the equilibrium population size of 
genotype 0 {K{r-m-s)lr vs. K{r-m)lr). 

In the second selection regime, we assume that this dif- 
ference in death rate depends on the population density in 
the surrounding grid cells. For bacteria located in a com- 
pletely empty neighborhood of / adjacent grid cells, allele 
0 entails a death rate that is changed by an amount Se 
(which can also take negative values). With increasing bac- 
terial density in this neighborhood, this selection coeffi- 
cient changes linearly until it reaches a value of Sj if the 
neighborhood is completely filled with bacteria. Hence, in- 
crease in death rate of genotype 0 is described by: 



Here Hi{i,j) is the set of / grid cells forming a neigh- 
borhood of grid cell (/,/); this also includes the grid cell 
(i,/) itself The term \Hi{i,J)\ denotes the size of Hi{i,j). 
The parameters used in the model and their units are 
listed in Table 1. 

We chose a realistic range of parameters for the pa- 
rameters that can be approximated experimentally, like 
the mutation and recombination rate. For other parame- 
ters, such as the rates of DNA diffusion and bacterial 
migration in a biofilm for which there are no exact ex- 
perimental estimates, the parameter space was screened. 
We assumed a mutation rate of 10"^. For the transform- 
ation rate, there is wide variation in experimentally ap- 
proximated transformation frequencies (10^^-10"^ for H. 
pylori [39] and 10"*^- 10"^ for Streptococcus pneumonia 
[40]). Here we used transformation rate of 10^^-10"^ 
(per bacterium per DNA fragment per unit time) in 
most simulations. We also studied the effect of higher 
transformation rates in some simulations. 

We examined the impact of transformation on adapta- 
tion in the course of colonization of a new environment. 
The initial bacterium had the genotype "0". We then 
followed the eco-evolutionary dynamics by numerically 
solving the above system of equations using the software 
package Mathematica v.8 (Wolfram Research Inc. [41]). 

Additional file 



Additional file 1: Figure SI. The effect of DNA decay on adaptation rate. 
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